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The next-nearest neighbor interaction (NNN) is included in a tight-binding calculation of the 
electronic spectrum and conductivity of doped graphene. As a result, we observe a wide variation 
of the conductivity behavior, since the Fermi energy and the resonance peak are not shifted by 
the same amount. Such effect can have a profound effect in the idea of explaining the minimal 
conductivity of graphene as a consequence of impurities or defects. Finally, we also estimate the 
mean free path and relaxation time due to resonant impurity scattering. 
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I. INTRODUCTION 

Research on graphene has seen a rush of publi- 
cations since its experimental discovery in 2004 
The vast fame of this carbon allotrope 0] is partly 
due to its amazing room-temperature transport 
properties as for example the high electronic 
[3] and thermal conductivity [5*1, profiling nano- 
devices based on graphene From the theoretical 
viewpoint, graphene is also astonishing because 
the charge carriers are described by massless Dirac 
fermions ^ [Tj] as a consequence of the crystal sym- 
metry. However, in the construction of electronic 
nano-devices, the use of pure graphene present 
some troubles. For example, the conductivity is 
difficult to manipulate by means of an external 
gate voltage, which is a desirable feature required 
to build a FET transistor. Such performance is 
related to the Klein paradox in relativistic quantum 
mechanics 7], or from a more standard outlook, 
as a consequence of the zero band gap. There are 
many proposals to solve the problem. For instance, 
by using quantum dots Q, a graphene nanomesh 
an external electromagnetic radiation source 
[TTI or by doping using impurities [11] • In fact, 
in a previous article we showed that impurities 
lead to a metal-insulator transition since a mobility 
edge appears near the Fermi energy [12|. Such 



prediction has been confirmed in graphene doped 
with H [l^, which opens the possibility to build 
graphene-based narrow gap semiconductors [T2j . 
Other groups have shown that graphene can present 
an n-typed semiconductor behavior when doped 
with N, Bi or Sb atoms; and a p-type semiconductor 
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behavior using B or Au atoms [flliBl- Still, there 
is much debate regarding the nature of the mobility 
transition, since in 2 dimensional (2D) scaling 
theory, it is predicted that all states are localized 
in the presence of a finite amount of disorder [l^ . 

The appearance of a mobility edge has its ori- 
gins in the presence of resonant states when few 
impurities are considered [Hi. Both type of states, 
localized and resonant, have an enhanced amplitude 
in the neighborhood of the impurity. Nevertheless, 
resonant states only trap electrons during a short 
time. Using a nearest-neighbor (NN) tight-binding 
model, resonant states have been reported near 
the Fermi energy [Ts! ]. Furthermore, an analytical 
approximate expression was found for the resonant 
energy as a function of the impurity energy using 
the Lifshitz equation [l^ [2^ . The main difference 
between such works is that Pereira et al. also 
noted that there is a slightly different between 
the resonance energy obtained from the Lifshitz 
equation, and the actual localization of the sharp 
resonance in the density of states (DOS) when the 
impurity energy is not so strong [l8| . It is necessary 
to remark that only str ong impurities are able to 
produce resonant states [21|. The presence of next- 
nearest neighbor interaction (NNN) shifts the Fermi 
energy and breaks the electron- hole symmetry ^22] . 
Including NNN interaction, the DOS displays a 
sharp peak when a vacancy is considered like an 
impurity in the lattice [23j . Moreover, this peak is 
smeared by the NNN. 

The central topic of this work is to emphasize 
the different behaviors in the electrical conductivity 
due to resonances when NNN interactions are 
included, as happens in carbon nanotubes [23 |. 
Usually, the NNN interaction is not taken into 



account in graphene tight-binding calculations 
[25I], so here we propose a systematic study of the 
subject. Such study is important because there is 
a debate concerning which mechanisms determine 
the charge carrier mobility [1^, [13], as well as the 
nature of the minimal conductivity 0]. As we will 
see, an impurity can produce a sharp peak or a 
smoothening effect in the electrical conductivity, 
depending on the charge doping, temperature, 
strength of the impurity scattering and the value of 
the NNN interaction. The interplay between such 
factors is subtle since for example, the Fermi level 
and the resonance energy are not shifted by the 
same amount when NNN interaction is included. It 
is worthwhile mentioning that the electrical conduc- 
tivity at high temperatures is determined basically 
by the electron-phonon interaction [2^, while here 
we discuss only scattering by impurities. Thus, 
our results are relevant only for low temperatures. 
However, this case is important to explain the weak 
temperature dependence of the conductivity, which 
is basically proportional to the carrier concentration 

i M- 

The layout of the work is the following. In sec- 
tion |TT] we describe the model and the perturbative 
approach used to calculate Green's function for a 
NNN tight binding Hamiltonian of doped graphene. 
Section Hill describes the calculation of Green's func- 
tion of pure graphene, which is used in Section IIVI 
to calculate the resonant energies. Section |V] con- 
tains the electrical conductivity calculations using 
the Kubo-Greenwood formula. Finally, in section 
I VII we present the conclusions. 



II. MODEL 

As a model, we consider a pure graphene tight- 
binding Hamiltonian with substitutional impurities 
at very low concentrations. Since there are no cor- 
relations between impurities and the impurity con- 
centration is very low, we can reduce the problem 
to a single localized impurity in a graphene lat- 
tice. The behavior for a given low concentration 
can be found by a simple implementation of the vir- 
tual crystal approximation (VGA) (30j. Also, we will 
use the fact that the graphene's honeycomb lattice 
is formed by two triangular interpenetrating sub- 
lattices, denoted A and B The corresponding 
tight-binding Hamiltonian is, 
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•Hi 



or Hi = e 



where ^ i'^t i) annihilates (creates) an electron 
with spin a (cr —'[, ],) on site i at position on the 
A sublattice (an equivalent definition is used for B 
sublattice), t (« 2.79eV) is the NN hopping energy, 
and t' (sa 0.68eV) is the NNN hopping energy {sil- 
e is the energy difference between a carbon atom 
and a foreign atom, and I is the impurity position. 

Usually, the resonances are characterized by look- 
ing at the Green's functions (G) of H. Expressing G 
as a perturbation series in terms of Gq (which is the 
Green's function corresponding to the unperturbed 
Hamiltonian Ho), a closed expression is obtained for 
the local density of states [LDOS] in the impurity 
site I [si, 



p{l-E) 



Po{l;E) 



\l-eGo{l,l;E)\^ 



(2) 



where Go{l,l;E) and po{l;E) are respectively the 
Green's function and LDOS on I site with Hq. 

The term |1 — eGo(^ l',E)\'^ cannot become zero for 
E within the band. However, under certain values 
of e, this term is near to zero for a given E E,-. 
Then, a sharp peak in the LDOS will emerge around 
Ej-. This E^ is associated with a resonant state inas- 
much as there is a different impurity energy level. 
If Im{Go(^ l',E)} is a slowly varying function of E 
(for E around i?i ), then the resonant energy will be 
given as a solution of the Lifshitz equation. 



l-eRe{Go(/,/;£;)}~0. 



(3) 



Furthermore, if the derivative of Re {Go(^, /; -E)} 
does not have a strong dependence of E near E^, 
then [H, 



1 



\l-eGo{lJ;E)\^ {E~E,Y+T^' 

where F corresponds to the width of the impurity 
resonance. 



H — Ho -l- Hi ; 



F = 



|Im{Go(/,/;^0}| 
|Re{G;(/, /;£;)} I 



(4) 
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Thus, the resonant state effect is sketched by its 
location, E^- in equation and its width, T in 
equation Those characteristics of the resonant 
state are inherited from the Green's function behav- 
ior, which is presented in the section below. 



e around E = 0, equation ([3]). This symmetry is 
broken when t' ^ 0, as seen in figure [T]('6j. This is 
due to the fact that real part of the Green's func- 
tion no longer cross the zero at E^ = 3^', this energy 
value matches with the zero value of the imaginary 
part. 



III. GREEN'S FUNCTION OF PURE 
GRAPHENE WITH NNN INTERACTION 

To solve the Lifshitz equation we need to ob- 
tain the Green's function for graphene with NNN 
interaction. Notice that analytical expressions are 
available only for the NN interaction [S^l, and not 
for the NNN interaction. The Green's function can 
be obtained from, 



Go(i?) 



- E 



1 



N ^ E + is-E(k) 



(5) 



where s <C 1 and £'(k) is the dispersion relationship 
of equation ([1]), and given by j23 |. 



i?±(k) = ±V3 + /(k)-i7(k) 



(6) 
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FIG. 1. (a) Green's functions with nearest neighbor in- 
teraction and (6) with next-nearest neighbor interaction. 
[The square mesh in the first Brillouin zone that was used 
to calculate Go{l, I', E) using the equation ([5]) is uniform 
and contains A*' = 7.5 x 10^ points, and s = 2 x 10~^.] 



/(k) = 2cos {y^ky 



-f 4 cos I —^kya 
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where the minus sign applies to the valence and the 
plus sign the conduction band. Around the Dirac 
point (K or K'), the momentum can be written as 
k = K -|- q, where q is a small vector. Then, equa- 
tion ^ up to second order is given by [11] , 



i?±(q) «3t'±^|q| 



± sm 
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(7) 



where a is the carbon-carbon distance (a « 1.42 A). 

To compute equation ([5]) we used a square mesh 
in the first Brillouin zone of the reciprocal space 
to evaluate the sum. The results are presented in 
figure [T] In figure we show the result when 
the NNN behavior is absent {t' = 0). The result is 
in excellent agreement with the analytical formula 
[3^ . Notice how at zero energy (corresponding to 
the Fermi energy for pure graphene, i?p), the imag- 
inary and real part of the Green's function cross at 
zero energy, resulting in a symmetrical behavior for 



IV. GREEN'S FUNCTION OF A SINGLE 
IMPURITY IN GRAPHENE 

Once the Green's function Go for pure graphene 
is known, we can compute G in order to describe 
the resonant states. In the following subsections 
we present the corresponding results. It is useful 
keeping in mind that the impurity effect is mainly 
weighted by the real part of Gq. 



A. Local Density of States (LDOS) 

The LDOS can be calculated using equation 
since 

Poil;E) = --lm{Goil,l;E)} . 

TT 

Using the calculated Go , and considering strong im- 
purities, i.e. e/t > 3, in figure [5] we can see that 
the LDOS exhibits a peak at certain resonant ener- 
gies for two combinations of impurities self-energies e 
and different NNN interaction t' . An evident char- 
acteristic is that the sharp location, i.e. the reso- 
nant energy, has a shift depending on the t' parame- 
ter, this is a consequence of the shift in the ordinate 
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axis of Re {Go{l, I; E)}, as observed in figured] An- 
other characteristic that corresponds to the sharpest 
LDOS behavior emerges when the E,- is near the 
Ep ■ From figure [U it is clear that the NNN interac- 
tion radicahy changes the resonance properties when 
compared with the NN case. Therefore, the NN in- 
teraction is not enough to describe the behavior of 
the doped system. To see this in more detail, let us 
calculate the position of the resonant energy and the 
resonance width. 
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FIG. 2. LDOS calculated using equation ((2]) for different 
values of e. andi'. [The parameters used to generate the 
graphs are the same as those used in Figure [T]] 



B. Resonant energy {Ei) 

In order to obtain E-^, we need to solve the Lif- 
shitz equation ([3]). As a result, in the figure [3] we 
present the energies £^r that satisfy equation ^ as 
a function of e for different sets of t' . The main ef- 
fect of the NNN interaction, is a shift proportional 
to t', as expected from the first correction to the NN 
interaction in equation (jT)). However, also the curva- 
ture of £'r(£) in figure[3]exhibits a slight difference as 
the NNN hopping energy varies. Additionally, in the 
same figure, we plot the Fermi energy E^ for pure 
graphene including NNN interaction for certain e at 
which the Fermi energy lies exactly at the resonance 
LDOS peak. In other words, the circles in figure |3] 
are the values of the parameters e and t' where the 
electronic properties are most affected, since elec- 
trons at the Fermi level have the exact energy of a 
resonant state and can thus be easily trapped for a 
certain amount of time around the impurity. 

According with Skrypnyk fio'], for the NN inter- 
action, the resonant states that are near E'^, in the 
asymptotic limit e — > ±cxd are given by. 
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FIG. 3. The = Er{£) curve defined by the Lifshitz 
equation ([Sjl, for different values of NNN interaction t' . 
The Fermi energy for pure graphene including NNN in- 
teraction, i?p, is identified by a circle. [The parameters 
used to generate the graphs are the same as those used 
in Figure [1]] 



In such limit, we can extend the previous expression 
to include the NNN interaction as follows, 
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(Notice that the numerical factor in the above ex- 
pression was not reported by Skrypnykv et al. [Tgj 
since their expression used for the dispersion rela- 
tion was not normalized). 

Although the previous expressions follows form a 
rigid translation of the spectrum with t' , the more 
realistic cases are those of small e, in which the ef- 
fects of t' are important, since the resonance peak is 
not shifted by the same amount. As we will see, this 
effect is important when one considers the effects on 
the electronic conductivity. 



C. Resonant width 

In figure 0] we present the influence of the NNN 
interaction t' on the resonance width F as a function 
of e^^. It is evident the influence of the NNN due 
to the asymmetry in the curves for t' ^ 0. Observe 
that the asymmetry leads to a reduced resonance 
width for e < when the NNN interaction increases. 
This means that the peak is shaper and thus, the 
lifetime of the resonance is increased. 

On the contrary, the opposite is observed for 
£ > 0; i.e., the lifetime of the electron near the 
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FIG. 4. Resonant width, F, given by equation Q as a 
function of the impurity effect, . The figure shows 
asymmetric curves due to the NNN interaction. [The 
parameters used to generate the graphs are the same as 
those used in Figure [T]] 



impurity is decreased by the NNN interaction. 

As was previously mentioned, the most notorious 
peak in the LDOS is located near the Fermi energy. 
The corresponding e value to that energy is clearly 
observed in figure SI which is the graph correspond- 
ing to equation (U). 



V. DC CONDUCTIVITY 

Resonances are important to the transport prop- 
erties, as for example, in the electrical or thermal 
conductivity. In this section, we evaluate the electri- 
cal conductivity {axx) taking into account resonant 
states and the NNN interaction. To do so, we use 
the Kubo-Greenwood formula expressed as [soj . 



^fi N 



(8) 



E+fi 



where, 



T{E)=TT{pxlm{G{E)}p.xlm{G{E)}} , 

and r^o ~ 3a^ is the primitive cell area, / is the 
Fermi-Dirac distribution and fj, is the chemical po- 
tential, which can be tuned by the external field (for 
example, with a voltage applied in the lattice), px 
is the momentum operator, given by the following 
commutator, 



Px = -^i^'2;J 



It is necessary to remark that here, the Hamiltonian 
operator includes the next-nearest neighbor interac- 
tion. Therefore, p^ inherits this interaction. After- 
math, Px can be written in terms of the momentum 
operator associated to the NN interaction as follows. 
Consider first the momentum px for NN, 



Px 



N 



7n Jx ' y ^ 



(9) 



I— I mGNN 



where we introduced the connectivity matrix defined 
as, 



W(to, n) = 



1 if m and n are NN 
otherwise 



Using this connectivity matrix, the Hamiltonian 
without perturbation including the NNN interaction 
can be rewritten as. 



Hq 



-tW -t'{W^ -31) , 



where I is the identity matrix. Taking the previous 
expression and using ^ , we obtain the correspond- 
ing operator for the NNN case. 



Px 



-tW-t'{W^ -31), x] 



j[pT^ + WpT 



(10) 



The previous expressions are effortlessly written 
in a computer program, in which we consider a 
low concentration of impurities, C, introduced 
adding the perturbation Hi with the same e at 
different sites taken at random with a uniform 
distribution. In figure [5] we show the conductivity 
calculated using the Kubo-Greenwood formula 
([8]). Considering different values of t' and as a 
function of the charge doping. Figure [5] was made 
at a fixed representative temperature, in this case 
/cbT = 0.025 eV, to highlight the main effects of 
the NNN interaction. Clearly, figure [S] exhibits the 
radical difference in the conductivity behavior due 
to the NNN interaction, since an smearing effect, as 
seen in figure [Slf'aJ, can appear, or as in figure [Slf^fc^, 
a sharp peak can be observed . In fact, if we look 
at the temperature behavior, we can also get very 
different behaviors of axx- These changes are due 
to a subtle interplay between the chemical doping, 
temperature, the NNN interaction and the impurity 
type. 



5 



t' = 0.00 
t' = 0.10 
i' = 0.2437 




e/t = 10 " 
-0.4 -0.2 0.2 0.4 

(fl - EF)/t 

(a) 



FIG. 5. Electrical conductivity calculated using the 
Kubo-Greenwood formula at low concentration of im- 
purities, C = 0.01. (a) Corresponds to e/t = 10 and (b) 
corresponds to e/t = —10. [All lattices have A'' ~ 10*.] 



To understand the diverse behaviors of (Txx, in 
figure [6] we show a sketch of the " building blocks" 
that appear in the Kubo-Greenwood, and how such 
blocks are modified by the considered parameters. 
On each panel of the graph, at the left we show 
the shape of the term df /dE, which corresponds 
to the "thermal selector". At T = 0, it becomes 
a delta function centered at the chemical potential 
(//), while at r 7^ has a width of the order of k^T. 
The position of df/dE on the energy axis can be 
externally modified by doping with charge carriers, 
resulting in different positions of the Fermi energy 
{Ep — fl) when compared with the equilibrium value 
of such energy, denoted by Ep . The second building 
block is T{E), which can be thought as a quantum 
selector of the transport channels. Since the Green's 
function of doped graphene can be written as. 



G ~ Gi 



Go\l)e{l\Go 



sGoiU) 



(11) 



(where the sum is carried over impurity sites) , T{E) 
has two behaviors. Near the resonant energies, G 
is dominated by the second term in Eq. (|11|) . and 
Go can be neglected. As a consequence, a peak ap- 
pears in the conductivity at the resonant energy, as 
shown in figure [H Far from the resonance, G ~ Gq. 
Then we recover the transmittance of pure graphene. 

Now we can study how the two building blocks 
interact to produce many different behaviors. First, 
it is clear that variations in fi and T can produce 
a peak in the conductivity if the thermal selector 
coincides with the resonance peak. So for example, 
the conductivity can be enhanced if for example, at 
a certain temperature the thermal selector begins 
to have an overlap over the peak of T{E). The 
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FIG. 6. Schematic diagram of the Kubo-Greenwood for- 
mula which explains the effect of the NNN interaction. 
The behavior of the two building blocks, the thermal se- 
lector of states df/dE and the trace T{E), are shown 
on each panel. The position of the resonance energy 
and Fermi energy are also shown. The following cases 
are considered: (a) e/t = 10 and (b) e/t — —10. Both 
graphics assume the case t' /t = 0.1. 



effect of the NNN interaction is very subtle since 
in principle one can expect a simple translation in 
energy of the spectrum. However, as it was said 
in the previous sections, the rigid translation of 
the spectrum is only valid at high e. According 
to our results, for realistic impurities there are 
deviations from such behavior, and thus, Ep and 
E-r are not rigidly translated, i.e., the distance 
\E-p — Er\ depends on the NNN interaction. Since 
the conductivity depends a lot on such factor due 
to the position of the thermal selector, the resulting 
effect of the NNN happens to be very important. 
One can see that in fact, the most important factor 
is the position of the resonant peak, which can be 
changed due to different kinds of impurities. Thus, 
one can expect a wide variability of the conduc- 
tivity due to resonant scattering, as is observed in 
nanotubes, wheres even similar samples can present 
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wide variations of the conductivity [2J| . 



Finally, if we consider the scattering term in Eq. 
(|lip , the impurities allow to define a relaxation time, 
Ts, as follows. 

where S is the scattering probability per unit time. 
S can be calculated by summing over all transition 
rates from \i) — ^ |/) , 

s.f = m - ff)w.f , 

using the Fermi golden rule, 

W,f^—\{fm)\^SiEf-E,), 

where 

1 r2 



After steps similar to those used to get the Kubo- 
Greenwood formula ISOll. we obtain, 



S 



2TTkBT 



(12) 



where. 



QiE) 



p\E) 



plil-Er){E-E,.y+r^ ■ 
Far from the resonance peak, Ey, and near Ep, 

then the scattering term Q{E) is given by. 



Q{E) 



1 1 



|£;-3i'P 



3n^pl{l;Er){3t'-E,)^ t* 



(13) 



At r = 0, the resulting relaxation time obtained 
from a straightforward evaluation of Eq. ([T2|) using 
(fT51) and remembering that Ep = 3t' , we obtain 



4C kuT 



\Et, - 3t' 



3h pI{1;E,) {it'-E,) 



Thus, the mean free path (/) is / « v-pTg, which 



goes as Ep"^ . 



VI. CONCLUSIONS 



We have studied the effects in the spectrum 
and electronic conductivity of low concentrations 
of impurities in graphene when the next-nearest 
neighbor interaction is considered in a tight-binding 
approximation. Although the electronic spectrum is 
basically similar to the case of pure nearest neighbor 
interaction, the conductivity is much more affected 
since the Fermi level and the resonance peak are 
not shifted by the same amount, resulting in a 
wide variability of the conductivity, as happens 
with carbon nanotubes [13]. As a consequence, the 
present study shows that is difficult to explain the 
minimum electrical conductivity of graphene at the 
Dirac point by only considering impurity scattering, 
and thus the evanescent waves approach seems to 
be a possible explanation Q. Also, we obtained 
that the relaxation time in graphene goes like -Ep^, 
resulting in a large mean free path. 
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